Leveraging in Vitro Models for Clinically Relevant Rare CYP2D6 Variants in Pharmacogenomics.

Cytochrome P450 2D6 (CYP2D6) is responsible for the metabolism of up to 20% of small-molecule drugs and therefore, may impact the safety and efficacy of medicines in broad therapeutic areas. CYP2D6 is highly polymorphic, and the frequency of variants can differ across racial and ethnic populations, significantly affecting enzymatic function and drug metabolism. However, rare variants of CYP2D6 present a unique challenge for academia, industry, and regulatory agencies alike due to the lack of feasibility of characterizing their clinical relevance in clinical trials, particularly in variants that exhibit population-specific frequencies in racial and ethnic groups that are poorly represented in clinical trials. Despite significant advancement in pharmacogenomics, the substrate specificity and related clinical relevance of these CYP2D6 rare variants remain largely unclear, and further efforts are warranted to characterize the burden of these variants on adverse drug reactions and drug efficacy. Thus, cell-based in vitro systems can be used to inform substrate-specific effects and the overall relevance of a rare variant. Liver microsomes, cell-based expression systems, ex vivo primary samples, and purified variant protein have all been used with various substrates to potentially predict the clinical impact of new substrates. In this review, we identify rare variants of CYP2D6 that demonstrate differences across races in prevalence and thus are often unassessed in clinical trials. Accordingly, we examine current pharmacogenomic in vitro models used to analyze the functional impact of these rare variants in a substrate-specific manner. SIGNIFICANCE STATEMENT: Variants of CYP2D6 play a clinically relevant role in drug metabolism, leading to potential safety and efficacy concerns. Although the influence of prevalent variants is often well characterized, rare variants are traditionally not included in clinical trials. This review captures the clinical relevance of rare variants in CYP2D6 by highlighting in vitro models that analyze their impact on the metabolism of CYP2D6 substrates.

The Allied Health Expansion Program: Rethinking how to prepare a workforce to enable improved public health outcomes.

Improvements in global public health require universal health care supported by a health workforce with competencies appropriate for local population needs-the right capabilities, in the right place, and at the right time. Health inequities persist in Tasmania, and Australia more broadly, most notably for those people living in rural and remote areas. The article describes the curriculum design thinking approach being used to codesign and develop a connected system of education and training to target intergenerational change in the allied health (AH) workforce capacity in Tasmania, and beyond. A curriculum design thinking process is engaging AH participant groups (faculty, AH professionals, and leaders across health, education, aged and disability sectors) in a series of focus groups and workshops. The design process deals with four questions: What is? What if? What wows? and What works? It also involves Discover, Define, Develop and Deliver phases that continue to inform the development of the new suite of AH education programs. The British Design Council's Double Diamond model is used to organize and interpret stakeholder input. During the initial design thinking discover phase, stakeholders identified four overarching problems: rurality, workforce challenges, graduate skill set shortfalls, and clinical placements and supervision. These problems are described in terms of relevance to the contextual learning environment in which AH education innovation is occurring. The develop phase of design thinking continues to involve working collaboratively with stakeholders to codesign potential solutions. Solutions to date include AH advocacy, a transformative visionary curriculum, and an interprofessional community-based education model. In Tasmania, innovative educational innovations are catalyzing attention and investment in the effective preparation of AH professionals for practice to deliver improved public health outcomes. A suite of AH education that is deeply networked and engaged with Tasmanian communities is being developed to drive transformational public health outcomes. These programs are playing an important role in strengthening the supply of allied health professionals with the right capabilities for metropolitan, regional, rural, and remote Tasmania. They are situated in a broader AH education and training strategy that supports the ongoing development of the AH workforce to better meet the therapy needs of people in Tasmanian communities.

Considerations for Use of Pharmacodynamic Biomarkers to Support Biosimilar Development - (III) A Randomized Trial with Interferon Beta-1a Products.

The US Food and Drug Administration (FDA) has taken steps to bring efficiency to the development of biosimilars, including establishing guidance for the use of pharmacokinetic and pharmacodynamic (PD) similarity study data without a comparative clinical study with efficacy end point(s). To better understand the potential role for PD biomarkers in biosimilar development and inform best practices for biomarker selection and analysis, we conducted a randomized, double-blinded, placebo-controlled, single-dose, parallel-arm clinical study in healthy participants. Eighty-four healthy participants (n = 12 per dose arm) received either placebo or one of three doses of either interferon ß-1a (7.5-30 µg) or pegylated interferon ß-1a (31.25-125 µg) to evaluate the maximum change from baseline and the baseline-adjusted area under the effect curve for the biomarkers neopterin in serum and myxovirus resistance protein 1 in blood. Both PD biomarkers increased following product administration with clear separation from baseline (neopterin: 3.4-fold and 3.9-fold increase for interferon ß-1a and pegylated interferon ß-1a, respectively; myxovirus resistance protein 1: 19.0-fold and 47.2-fold increase for interferon ß-1a and pegylated interferon ß-1a, respectively). The dose-response curves support that therapeutic doses were adequately sensitive to detect differences in both PD biomarkers for consideration in a PD similarity study design. Because baseline levels of both biomarkers are low compared with on-treatment values, there was little difference in using PD measures adjusted to baseline compared with the results without baseline adjustment. This study illustrates potential methodologies for evaluating PD biomarkers and an approach to address information gaps when limited information is publicly available for one or more PD biomarkers.

Evaluating the Utility of Proteomics for the Identification of Circulating Pharmacodynamic Biomarkers of IFNß-1a Biologics.

Proteomics has the potential to identify pharmacodynamic (PD) biomarkers for similarity assessment of proposed biosimilars without relying on clinical efficacy end points. In this study, with 36 healthy participants randomized to therapeutic doses of interferon-beta 1a products (IFNß-1a) or pegylated-IFNß-1a (pegIFNß-1a) approved to treat multiple sclerosis or placebo, we evaluated the utility of a proteomic assay that profiles > 7,000 plasma proteins. IFNß-1a and pegIFNß-1a resulted in 248 and 528 differentially expressed protein analytes, respectively, between treatment and placebo groups over the time course. Thirty-one proteins were prioritized based on a maximal fold change ≥ 2 from baseline, baseline adjusted area under the effect curve (AUEC) and overlap between the 2 products. Of these, the majority had a significant AUEC compared with placebo in response to either product; 8 proteins showed > 4-fold maximal change from baseline. We identified previously reported candidates, beta-2microglobulin and interferon-induced GTP-binding protein (Mx1) with ~ 50% coefficient of variation (CV) for AUEC, and many new candidates (including I-TAC, C1QC, and IP-10) with CVs ranging from 26%-129%. Upstream regulator analysis of differentially expressed proteins predicted activation of IFNß1 signaling as well as other cytokine, enzyme, and transcription signaling networks by both products. Although independent replication is required to confirm present results, our study demonstrates the utility of proteomics for the identification of individual and composite candidate PD biomarkers that may be leveraged to support clinical pharmacology studies for biosimilar approvals, especially when biologics have complex mechanisms of action or do not have previously characterized PD biomarkers.

Expanding Approved Patient Populations for Rare Disease Treatment Using In Vitro Data.

In vitro cell-based data can be used to support the extension of pharmaceutical approval to patient subsets with unique genetic variants. A set of conditions should be satisfied to support the extension of approval. The disease mechanism should be well described, and the impact of variants on protein function should be reasonably understood. The incidence data should show that clinical trials for the variants in question are not practical. The overall safety and efficacy of the drug should be clear in adequate and well-controlled clinical trials. The clinical trial should include patients found to be responders and nonresponders so that both positive and negative predictive power of the in vitro assay may be measured. The mechanism of action of the drug should be clearly defined and should be consistent with the disease mechanism. The assay system should be qualified, including the following points: (i) each variant construct should be confirmed by bidirectional sequencing; (ii) the in vitro assay should directly measure the variant protein function in comparison with the reference protein; (iii) the assay should be formally validated to the extent possible, clearly demonstrating precision, reproducibility, and sensitivity used to support the efficacy claim; and (iv) the primary data should be available for inspection and analytical validation. The overall goal is a robust and validated cell-based system that can be shown to predict the outcome of targeted therapy.

New science, drug regulation, and emergent public health issues: The work of FDA's division of applied regulatory science.

The U.S. Food and Drug Administration (FDA) Division of Applied Regulatory Science (DARS) moves new science into the drug review process and addresses emergent regulatory and public health questions for the Agency. By forming interdisciplinary teams, DARS conducts mission-critical research to provide answers to scientific questions and solutions to regulatory challenges. Staffed by experts across the translational research spectrum, DARS forms synergies by pulling together scientists and experts from diverse backgrounds to collaborate in tackling some of the most complex challenges facing FDA. This includes (but is not limited to) assessing the systemic absorption of sunscreens, evaluating whether certain drugs can convert to carcinogens in people, studying drug interactions with opioids, optimizing opioid antagonist dosing in community settings, removing barriers to biosimilar and generic drug development, and advancing therapeutic development for rare diseases. FDA tasks DARS with wide ranging issues that encompass regulatory science; DARS, in turn, helps the Agency solve these challenges. The impact of DARS research is felt by patients, the pharmaceutical industry, and fellow regulators. This article reviews applied research projects and initiatives led by DARS and conducts a deeper dive into select examples illustrating the impactful work of the Division.

Genome-wide association study of café-au-lait macule number in neurofibromatosis type 1.

BACKGROUND: Neurofibromatosis type 1 (NF1) is a tumor-predisposition disorder that arises due to pathogenic variants in tumor suppressor NF1. NF1 has variable expressivity that may be due, at least in part, from heritable elements such as modifier genes; however, few genetic modifiers have been identified to date. METHODS: In this study, we performed a genome-wide association analysis of the number of café-au-lait macules (CALM) that are considered a tumor-like trait as a clinical phenotype modifying NF1. RESULTS: A borderline genome-wide significant association was identified in the discovery cohort (CALM1, N = 112) between CALM number and rs12190451 (and rs3799603, r2  = 1.0; p = 7.4 × 10-8 ) in the intronic region of RPS6KA2. Although, this association was not replicated in the second cohort (CALM2, N = 59) and a meta-analysis did not show significantly associated variants in this region, a significant corroboration score (0.72) was obtained for the RPS6KA2 signal in the discovery cohort (CALM1) using Complementary Pairs Stability Selection for Genome-Wide Association Studies (ComPaSS-GWAS) analysis, suggesting that the lack of replication may be due to heterogeneity of the cohorts rather than type I error. CONCLUSION: rs12190451 is located in a melanocyte-specific enhancer and may influence RPS6KA2 expression in melanocytes-warranting further functional studies.

Associations Between Prediagnostic Concentrations of Circulating Sex Steroid Hormones and Esophageal/Gastric Cardia Adenocarcinoma Among Men.

Background: Esophageal adenocarcinoma (EA) and gastric cardia adenocarcinoma (GCA) are characterized by a strong male predominance. Concentrations of sex steroid hormones have been hypothesized to explain this sex disparity. However, no prospective population-based study has examined sex steroid hormones in relation to EA/GCA risk. Thus, we investigated whether prediagnostic circulating sex steroid hormone concentrations were associated with EA/GCA in a nested case-control study drawn from participants in three prospective cohort studies. Methods: Using gas chromatography-mass spectrometry (GC-MS) and electrochemiluminescence immunoassay, we quantitated sex steroid hormones and sex hormone binding globulin, respectively, in serum from 259 EA/GCA male case participants and 259 matched male control participants from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial, Alpha-Tocopherol, Beta-Carotene Cancer Prevention Study, and Cancer Prevention Study II Nutrition Cohort. Multivariable conditional logistic regression was used to calculate odds ratios (ORs) and 95% confidence intervals (CIs) for associations between circulating hormones and EA/GCA risk. All statistical tests were two-sided. Results: Higher concentrations of dehydroepiandrosterone (DHEA) were associated with a 38% decreased risk of EA/GCA (OR per unit increase in log2 DHEA = 0.62, 95% CI = 0.47 to 0.82, Ptrend = .001). Higher estradiol concentrations were associated with a 34% reduced risk of EA/GCA (OR = 0.66, 95% CI = 0.45 to 0.98, Ptrend = .05), and the association with free estradiol was similar. No other associations between baseline hormone concentrations and future EA/GCA risk were observed. Conclusions: This study provides the first evidence that higher concentrations of circulating DHEA, estradiol, and free estradiol may be associated with lower risks of EA/GCA in men.

Agnostic Pathway/Gene Set Analysis of Genome-Wide Association Data Identifies Associations for Pancreatic Cancer.

BACKGROUND: Genome-wide association studies (GWAS) identify associations of individual single-nucleotide polymorphisms (SNPs) with cancer risk but usually only explain a fraction of the inherited variability. Pathway analysis of genetic variants is a powerful tool to identify networks of susceptibility genes. METHODS: We conducted a large agnostic pathway-based meta-analysis of GWAS data using the summary-based adaptive rank truncated product method to identify gene sets and pathways associated with pancreatic ductal adenocarcinoma (PDAC) in 9040 cases and 12 496 controls. We performed expression quantitative trait loci (eQTL) analysis and functional annotation of the top SNPs in genes contributing to the top associated pathways and gene sets. All statistical tests were two-sided. RESULTS: We identified 14 pathways and gene sets associated with PDAC at a false discovery rate of less than 0.05. After Bonferroni correction (P ≤ 1.3 × 10-5), the strongest associations were detected in five pathways and gene sets, including maturity-onset diabetes of the young, regulation of beta-cell development, role of epidermal growth factor (EGF) receptor transactivation by G protein-coupled receptors in cardiac hypertrophy pathways, and the Nikolsky breast cancer chr17q11-q21 amplicon and Pujana ATM Pearson correlation coefficient (PCC) network gene sets. We identified and validated rs876493 and three correlating SNPs (PGAP3) and rs3124737 (CASP7) from the Pujana ATM PCC gene set as eQTLs in two normal derived pancreas tissue datasets. CONCLUSION: Our agnostic pathway and gene set analysis integrated with functional annotation and eQTL analysis provides insight into genes and pathways that may be biologically relevant for risk of PDAC, including those not previously identified.

Association between circulating levels of sex steroid hormones and esophageal adenocarcinoma in the FINBAR Study.

BACKGROUND: Esophageal adenocarcinoma (EA) is characterized by a strong male predominance. Sex steroid hormones have been hypothesized to underlie this sex disparity, but no population-based study to date has examined this potential association. METHODS: Using mass spectrometry and ELISA, we quantitated sex steroid hormones and sex hormone binding globulin, respectively, in plasma from males- 172 EA cases and 185 controls-within the Factors Influencing the Barrett/Adenocarcinoma Relationship (FINBAR) Study, a case-control investigation conducted in Northern Ireland and Ireland. Multivariable adjusted logistic regression was used to calculate odds ratios (ORs) and 95% confidence intervals (CIs) for associations between circulating hormones and EA. RESULTS: Higher androgen:estrogen ratio metrics were associated with increased odds of EA (e.g., testosterone:estradiol ratio ORQ4 v. Q1 = 2.58, 95%CI = 1.23-5.43; Ptrend = 0.009). All estrogens and androgens were associated with significant decreased odds of EA. When restricted to individuals with minimal to no decrease in body mass index, the size of association for the androgen:estrogen ratio was not greatly altered. CONCLUSIONS: This first study of sex steroid hormones and EA provides tentative evidence that androgen:estrogen balance may be a factor related to EA. Replication of these findings in prospective studies is needed to enhance confidence in the causality of this effect.

Association of high-evidence gastric cancer susceptibility loci and somatic gene expression levels with survival.

Eleven high-evidence single-nucleotide polymorphisms (SNPs) at nine loci for gastric cancer (GC) risk were reported, but their associations with survival remain unknown. In this study, we examined associations between SNP and GC survival by anatomic location and histology among 1147 incident cases from the Shanxi Upper Gastrointestinal Genetics Project. We further examined whether SNPs were expression quantitative trait loci in normal and tumor gastric tissues, and whether tumor versus normal somatic mRNA differences in 126 cases were associated with survival. No SNPs were associated with GC survival overall. However, subtype-specific associations were observed for gastric cardia adenocarcinomas at MUC1/TRIM46/1q22 rs2070803 [HRAA versus GA+GG = 2.16; 95% confidence interval (CI) = 1.24-3.78; P = 0.0068] and LTA/TNF/6p21.33 rs1799724 (HRTT+CT versus CC = 1.30; 95% CI = 1.07-1.57; P = 0.0077), and for diffuse-type GC at PSCA/8q24.3 rs2294008 (HRTT versus CT+CC = 1.99; 95% CI = 1.33-2.97; P = 7.8E-04). Rs2294008T was a cis-expression quantitative trait loci for PSCA, upregulating mRNA in normal gastric (ß = 0.60; P = 5.7E-21) and GC (ß = 0.30; P = 0.0089) tissues. Cases in the highest quartile (the smallest downregulation of tumor PSCA) had shortest survival than cases with the most downregulated PSCA (median survival of 0.47 years in the highest quartile versus 3.73 years in the lowest quartile; hazard ratio = 9.70; 95% CI = 2.46-38.4; P = 0.0012). Less striking effects for mRNA levels were observed for MTX1 at 1q22 in gastric cardia adenocarcinoma and for JRK at 8q24.3 in diffuse GC. Our results suggest three high-evidence GC risk loci have prognostic importance in GC subtypes. Future studies in well-characterized independent populations are warranted to validate our findings and further investigate the clinical utility of these variants in predicting GC prognosis.

GWAS follow-up study of esophageal squamous cell carcinoma identifies potential genetic loci associated with family history of upper gastrointestinal cancer.

Based on our initial genome-wide association study (GWAS) on esophageal squamous cell carcinoma (ESCC) in Han Chinese, we conducted a follow-up study to examine the single nucleotide polymorphisms (SNPs) associated with family history (FH) of upper gastrointestinal cancer (UGI) cancer in cases with ESCC. We evaluated the association between SNPs and FH of UGI cancer among ESCC cases in a stage-1 case-only analysis of the National Cancer Institute (NCI, 541 cases with FH and 1399 without FH) and Henan GWAS (493 cases with FH and 869 without FH) data (discovery phase). The top SNPs (or their surrogates) from discovery were advanced to a stage-2 evaluation in additional Henan subjects (2801 cases with FH and 3136 without FH, replication phase). A total of 19 SNPs were associated with FH of UGI cancer in ESCC cases with P < 10-5 in the stage-1 meta-analysis of NCI and Henan GWAS data. In stage-2, the association for rs79747906 (located at 18p11.31, P = 5.79 × 10-6 in discovery) was replicated (P = 0.006), with a pooled-OR of 1.59 (95%CI: 1.11-2.28). We identified potential genetic variants associated with FH of UGI cancer. Our findings may provide important insights into new low-penetrance susceptibility regions involved in the susceptibility of families with multiple UGI cancer cases.

Inherited variation in circadian rhythm genes and risks of prostate cancer and three other cancer sites in combined cancer consortia.

Circadian disruption has been linked to carcinogenesis in animal models, but the evidence in humans is inconclusive. Genetic variation in circadian rhythm genes provides a tool to investigate such associations. We examined associations of genetic variation in nine core circadian rhythm genes and six melatonin pathway genes with risk of colorectal, lung, ovarian and prostate cancers using data from the Genetic Associations and Mechanisms in Oncology (GAME-ON) network. The major results for prostate cancer were replicated in the Prostate, Lung, Colorectal and Ovarian (PLCO) cancer screening trial, and for colorectal cancer in the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). The total number of cancer cases and controls was 15,838/18,159 for colorectal, 14,818/14,227 for prostate, 12,537/17,285 for lung and 4,369/9,123 for ovary. For each cancer site, we conducted gene-based and pathway-based analyses by applying the summary-based Adaptive Rank Truncated Product method (sARTP) on the summary association statistics for each SNP within the candidate gene regions. Aggregate genetic variation in circadian rhythm and melatonin pathways were significantly associated with the risk of prostate cancer in data combining GAME-ON and PLCO, after Bonferroni correction (ppathway < 0.00625). The two most significant genes were NPAS2 (pgene = 0.0062) and AANAT (pgene = 0.00078); the latter being significant after Bonferroni correction. For colorectal cancer, we observed a suggestive association with the circadian rhythm pathway in GAME-ON (ppathway = 0.021); this association was not confirmed in GECCO (ppathway = 0.76) or the combined data (ppathway = 0.17). No significant association was observed for ovarian and lung cancer. These findings support a potential role for circadian rhythm and melatonin pathways in prostate carcinogenesis. Further functional studies are needed to better understand the underlying biologic mechanisms.

Rare Germline Copy Number Variations and Disease Susceptibility in Familial Melanoma.

Mounting evidence suggests that copy number variations (CNVs) can contribute to cancer susceptibility. The main goal of this study was to evaluate the role of germline CNVs in melanoma predisposition in high-risk melanoma families. We used genome-wide tiling comparative genomic hybridization and single nucleotide polymorphism arrays to characterize CNVs in 335 individuals (240 melanoma cases) from American melanoma-prone families (22 with germline CDKN2A or CDK4 mutations). We found that the global burden of overall CNVs (or deletions or duplications separately) was not significantly associated with case-control or CDKN2A/CDK4 mutation status after accounting for the familial dependence. However, we identified several rare CNVs that either involved known melanoma genes (e.g., PARP1, CDKN2A) or cosegregated with melanoma (duplication on 10q23.23, 3p12.2 and deletions on 8q424.3, 2q22.1) in families without mutations in known melanoma high-risk genes. Some of these CNVs were correlated with expression changes in disrupted genes based on RNASeq data from a subset of melanoma cases included in the CNV study. These results suggest that rare cosegregating CNVs may influence melanoma susceptibility in some melanoma-prone families and genes found in our study warrant further evaluation in future genetic analyses of melanoma.

A Powerful Procedure for Pathway-Based Meta-analysis Using Summary Statistics Identifies 43 Pathways Associated with Type II Diabetes in European Populations.

Meta-analysis of multiple genome-wide association studies (GWAS) has become an effective approach for detecting single nucleotide polymorphism (SNP) associations with complex traits. However, it is difficult to integrate the readily accessible SNP-level summary statistics from a meta-analysis into more powerful multi-marker testing procedures, which generally require individual-level genetic data. We developed a general procedure called Summary based Adaptive Rank Truncated Product (sARTP) for conducting gene and pathway meta-analysis that uses only SNP-level summary statistics in combination with genotype correlation estimated from a panel of individual-level genetic data. We demonstrated the validity and power advantage of sARTP through empirical and simulated data. We conducted a comprehensive pathway-based meta-analysis with sARTP on type 2 diabetes (T2D) by integrating SNP-level summary statistics from two large studies consisting of 19,809 T2D cases and 111,181 controls with European ancestry. Among 4,713 candidate pathways from which genes in neighborhoods of 170 GWAS established T2D loci were excluded, we detected 43 T2D globally significant pathways (with Bonferroni corrected p-values < 0.05), which included the insulin signaling pathway and T2D pathway defined by KEGG, as well as the pathways defined according to specific gene expression patterns on pancreatic adenocarcinoma, hepatocellular carcinoma, and bladder carcinoma. Using summary data from 8 eastern Asian T2D GWAS with 6,952 cases and 11,865 controls, we showed 7 out of the 43 pathways identified in European populations remained to be significant in eastern Asians at the false discovery rate of 0.1. We created an R package and a web-based tool for sARTP with the capability to analyze pathways with thousands of genes and tens of thousands of SNPs.

Whole exome sequencing in families at high risk for Hodgkin lymphoma: identification of a predisposing mutation in the KDR gene.

Hodgkin lymphoma shows strong familial aggregation but no major susceptibility genes have been identified to date. The goal of this study was to identify high-penetrance variants using whole exome sequencing in 17 Hodgkin lymphoma prone families with three or more affected cases or obligate carriers (69 individuals), followed by targeted sequencing in an additional 48 smaller HL families (80 individuals). Alignment and variant calling were performed using standard methods. Dominantly segregating, rare, coding or potentially functional variants were further prioritized based on predicted deleteriousness, conservation, and potential importance in lymphoid malignancy pathways. We selected 23 genes for targeted sequencing. Only the p.A1065T variant in KDR (kinase insert domain receptor) also known as VEGFR2 (vascular endothelial growth factor receptor 2) was replicated in two independent Hodgkin lymphoma families. KDR is a type III receptor tyrosine kinase, the main mediator of vascular endothelial growth factor induced proliferation, survival, and migration. Its activity is associated with several diseases including lymphoma. Functional experiments have shown that p.A1065T, located in the activation loop, can promote constitutive autophosphorylation on tyrosine in the absence of vascular endothelial growth factor and that the kinase activity was abrogated after exposure to kinase inhibitors. A few other promising mutations were identified but appear to be "private". In conclusion, in the largest sequenced cohort of Hodgkin lymphoma families to date, we identified a causal mutation in the KDR gene. While independent validation is needed, this mutation may increase downstream tumor cell proliferation activity and might be a candidate for targeted therapy.

MEGSA: A Powerful and Flexible Framework for Analyzing Mutual Exclusivity of Tumor Mutations.

The central challenges in tumor sequencing studies is to identify driver genes and pathways, investigate their functional relationships, and nominate drug targets. The efficiency of these analyses, particularly for infrequently mutated genes, is compromised when subjects carry different combinations of driver mutations. Mutual exclusivity analysis helps address these challenges. To identify mutually exclusive gene sets (MEGS), we developed a powerful and flexible analytic framework based on a likelihood ratio test and a model selection procedure. Extensive simulations demonstrated that our method outperformed existing methods for both statistical power and the capability of identifying the exact MEGS, particularly for highly imbalanced MEGS. Our method can be used for de novo discovery, for pathway-guided searches, or for expanding established small MEGS. We applied our method to the whole-exome sequencing data for 13 cancer types from The Cancer Genome Atlas (TCGA). We identified multiple previously unreported non-pairwise MEGS in multiple cancer types. For acute myeloid leukemia, we identified a MEGS with five genes (FLT3, IDH2, NRAS, KIT, and TP53) and a MEGS (NPM1, TP53, and RUNX1) whose mutation status was strongly associated with survival (p = 6.7 × 10(-4)). For breast cancer, we identified a significant MEGS consisting of TP53 and four infrequently mutated genes (ARID1A, AKT1, MED23, and TBL1XR1), providing support for their role as cancer drivers.

Pathway, in silico and tissue-specific expression quantitative analyses of oesophageal squamous cell carcinoma genome-wide association studies data.

BACKGROUND: Oesophageal cancer is the fourth leading cause of cancer death in China where essentially all cases are histologically oesophageal squamous cell carcinoma (ESCC). Agnostic pathway-based analyses of genome-wide association study (GWAS) data combined with tissue-specific expression quantitative trait loci (eQTL) analysis and publicly available functional data can identify biological pathways and/or genes enriched with functionally-relevant disease-associated variants. METHOD: We used the adaptive multilocus joint test to analyse 1827 pathways containing 6060 genes using GWAS data from 1942 ESCC cases and 2111 controls with Chinese ancestry. We examined the function of risk alleles using in silico and eQTL analyses in oesophageal tissues. RESULTS: Associations with ESCC risk were observed for 36 pathways predominantly involved in apoptosis, cell cycle regulation and DNA repair and containing known GWAS-associated genes. After excluding genes with previous GWAS signals, candidate pathways (and genes) for ESCC risk included taste transduction (KEGG_hsa04742; TAS2R13, TAS2R42, TAS2R14, TAS2R46,TAS2R50), long-patch base excision repair (Reactome_pid; POLD2) and the metabolics pathway (KEGG_hsa01100; MTAP, GAPDH, DCTD, POLD2, AMDHD1). We identified and validated CASP8 rs13016963 and IDH2 rs11630814 as eQTLs, and CASP8 rs3769823 and IDH2 rs4561444 as the potential functional variants in high-linkage disequilibrium with these single nucleotide polymorphisms (SNPs), respectively. Further, IDH2 mRNA levels were down-regulated in ESCC (tumour:normal-fold change = 0.69, P = .75E-14). CONCLUSION: Agnostic pathway-based analyses and integration of multiple types of functional data provide new evidence for the contribution of genes in taste transduction and metabolism to ESCC susceptibility, and for the functionality of both established and new ESCC risk-related SNPs.